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By selectively probing the center of a trapped gas, we measure the local, or homogeneous, contact 
of a unitary Fermi gas as a function of temperature. Tan's contact, C, is proportional to the deriva- 
tive of the energy with respect to the interaction strength, and is thus an essential thermodynamic 
quantity for a gas with short-range correlations. Theoretical predictions for the temperature depen- 
dence of C differ substantially, especially near the superfluid transition, T c , where C is predicted 
to either sharply decrease, sharply increase, or change very little. For T/Tf > 0.4, our measure- 
ments of the homogeneous gas contact show a gradual decrease of C with increasing temperature, 
as predicted by theory. We observe a sharp decrease in C at T/Tf = 0.16, which may be due to the 
superfluid phase transition. While a sharp decrease in C below T c is predicted by some many-body 
theories, we find that none of the predictions fully accounts for the data. 
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The collective behavior of an ensemble of strongly in- 
teracting fermions is central to many physical systems 
including liquid 3 He, high-T c superconductors, quark- 
gluon plasma, neutron stars, and ultracold Fermi gases. 
However, theoretical understanding of strongly interact- 
ing fermions is challenging due to the many-body nature 
of the problem and the fact that there is no obvious small 
parameter for a perturbative analysis. Therefore, in or- 
der to establish the validity and applicability of theoret- 
ical approaches, it is essential to compare them against 
experimental results. Ultracold atomic Fermi gases are 
ideal for this purpose, as they provide excellent control- 
lability, reproducibility, and unique detection methods 
[U [3] ■ In particular, changing the magnetic field in the 
vicinity of a Feshbach resonance enables precise control 
of the interactions, which are characterized by the s-wave 
scattering length 3J. On resonance, the scattering length 
diverges and the behavior of the unitary gas no longer de- 
pends on it. Testing theories in this regime is especially 
desirable. 

An outstanding issue for the unitary Fermi gas is the 
nature of the normal state just above the transition tem- 
perature, T c , for a superfluid of paired fermions. Some 
theories of strongly interacting Fermi gases (BCS-BEC 
crossover theories) predict that the normal state is not 
the ubiquitous Fermi liquid but instead involves incoher- 
ent fermion pairing (preformed pairs) in what has been 
termed the pseudogap state @]. It has been suggested 
that the pseudogap state affects the temperature depen- 
dence of a quantity called Tan's contact [ 5 j. The con- 
tact, which is a measure of the short-range correlation 
function, has been shown to be an essential thermody- 
namic parameter for ensembles with short-range interac- 
tions [fJHTTj . The contact connects many seemingly unre- 
lated quantities through a set of universal relations that 
are valid for any temperature, any interaction strength, 
and any phase of the system. While the value of the 
contact, as well as many of these relations, were tested 



successfully at low temperature [j"2"HT5] , there are signifi- 
cant discrepancies among theories on how the contact of 
a unitary homogeneous Fermi gas depends on tempera- 
ture, especially around T c USHTrj] . The temperature 
dependence of the contact was recently measured for a 
trapped unitary Fermi gas |20j . However, for the trapped 
gas, averaging over the inhomogeneous density distribu- 
tion washes out any temperature-dependent features, and 
the measurement was unable to differentiate between the- 
oretical models. Here we present a measurement of the 
homogeneous contact, which can be directly compared to 
the predictions of different many-body theories. 

We perform the experiments with an optically trapped 
ultracold gas of 40 K atoms in an equal mixture of the 
\f,m f ) = |9/2, —9/2) and |9/2, —7/2) spin states, where 
/ is the quantum number denoting the total atomic 
spin and m/ is its projection |21j . We determine the 
contact by combining rf spectroscopy with a recently 
demonstrated technique to probe the local properties of a 
trapped gas [22]. Probing atoms locally is accomplished 
by intersecting two perpendicularly propagating hollow 
light beams that optically pump atoms at the edge of 
the cloud into a spin state that is dark to the detection 
(sec figure [l^,). The experimental sequence is depicted in 
figure [T]d. The magnetic field is ramped adiabatically to 
the Feshbach resonance and kept at this value for 2 ms 
before abruptly shutting off the trapping potential. Be- 
fore the potential is shut off, the hollow beams are pulsed 
on, followed by the rf pulse, which transfers a small frac- 
tion of the atoms in the occupied |9/2, —7/2) state to the 
initially unoccupied |9/2, —5/2) state (which is weakly in- 
teracting with the other two spin states) . We detect these 
atoms using absorption imaging after 3 ms of expansion. 
The temperature of the gas is varied by changing the fi- 
nal depth of the optical dipole trap in the evaporation 
process [21]. The number of atoms per spin state af- 
ter the evaporation ranges from 50,000 to 220,000. For 
the data presented in this paper, the radial trapping fre- 
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FIG. 1. Schematics of the experiment, (a) We probe the 
center of the gas by optically pumping the outer parts of 
the cigar-shaped cloud to a dark state using two intersect- 
ing second-order Laguerre-Gaussian beams [22] . By changing 
the beams' power, we control the fraction of atoms probed. 
The power of the two beams is set such that the number of 
atoms optically pumped by each beam is about the same, (b) 
The magnetic field is ramped from 203.4 G, where the atoms 
are initially prepared, to the Feshbach resonance. The hollow 
light beams are turned on 280 fis before trap release; initially, 
the beam that propagates perpendicular to the long axis of 
the cloud is pulsed on for 10 [is followed by 40 fis of the sec- 
ond beam. The line shape is measured using an rf pulse with 
a total duration of 100 /is and a gaussian field envelope with 
a = 17 fis, centered 180 /^s before trap release. The cloud ex- 
pands for 3 ms before being detected by absorption imaging. 
To improve the signal-to-noise ratio, we remove the remaining 
atoms from the |9/2, —9/2) and j 9/2, -7/2) states and then 
transfer the outcoupled atoms in the 1 9/2 , —5/2) state to the 
1 9/2, — 9/2) state, where we image on the cycling transition 
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FIG. 2. An rf line shape for the unitary Fermi gas at T/Tf = 
0.25 with 30% of the atoms probed. The solid (red) line is a 
fit to Eq. dll with the normalization Y(v)dv = 1/2, due 
to the 50% — 50% spin mixture. The inset shows the same 
data multiplied by 2 3 ^ 2 ty 2 u 3 ^ 2 . We make sure the rf pulse 
induces only a small perturbation, by setting its power to 
well below the value where we see the onset of saturation of 
the number of outcoupled atoms [24]. The measurement at 
different frequencies is done with different rf powers, and when 
analyzing the data, we linearly scale the measured number 
of atoms outcoupled at each frequency to correspond to a 
common rf power. 



quency, uj t , ranges from 2tt x 200 Hz to 2tt x 410 Hz, while 
the axial trapping frequency, u) z , ranges from 2-7T x 19 Hz 
to 2vr x 25 Hz. 

The contact is extracted from a measurement of the 
rf line shape T(v) [12], where T(v) is the rate of atoms 
transferred from one of the two interacting spin states 
to a third state, by an rf pulse centered at a frequency 
detuning v. A representative data set, where the hollow 
light beams were used to select the central 30% of the 
atom cloud, is shown in figure [2j For each line shape, 
we take data at 30 different detunings between —16 kHz 
and +116 kHz, where v = is defined as the single- 
particle transition frequency between the 1 9/2, —7/2) and 
1 9/2, — 5/2) states (measured for a spin polarized gas in 
the 1 9/2, — 7/2) state). The highest frequency typically 
corresponds to approximately 13 Ep/h. The high fre- 
quency tail of the rf line shape is predicted to scale as 
j/~ 3 / 2 , with the contact connecting the amplitude of the 
high frequency tail through (see Ref. [TT| and references 
therein) : 



I» 



= C/(Nk F ) 



for 



(1) 



where N is the total number of atoms, and Hkp is the 
Fermi momentum, and v is the rf detuning in units of the 
Fermi energy, Ep/h, with h being the Planck constant 
(2irh = h). The inset of figure [2] shows T(v) multiplied 



by 2 3 / 2 7t 2 j/ 3 / 2 , where we observe a plateau for frequen- 
cies higher than 5 Ep/h. We extract the contact by fit- 
ting the measured T(i>) for v > 5 Ep/h to Eq.Q (solid 
line in figure |2j). For the normalization, we integrate the 
line shape, including the tail, up to v — h/mr^ s , where 
r c ff is the effective range of the interaction [3] (which is 
approximately 300Ep/h). 

The main result of the paper, namely the homoge- 
neous contact versus the temperature, is presented in 
figure [3] The contact is normalized to the average kp 
of the probed sample, and temperature is given in terms 
of T/Tp, with Tp being the average Fermi temperature 
of the probed sample (we explain later how we determine 
these quantities). The data shows a monotonic decrease 
of the contact with increasing temperature from a maxi- 
mum value of 3.3 Nk F . For T /Tp = 0.16, at the edge of 
our experimentally attainable temperatures, we observe 
a sharp decrease of the contact to about 2.6 Nkp. In 
figure [3 we compare our data with several theoretical 
models [18] and a quantum Monte-Carlo (QMC) simula- 
tion |19] . The many-body theories are in the framework 
of the t-matrix approximation, differing by their choice of 
the diagrammatic expansion, the particle-particle prop- 
agator, and the self-energy. For T/Tp > 0.4, the differ- 
ences between the theoretical models are small, and the 
predictions all lie within the uncertainty of the data. As 
expected, at higher temperatures (T/Tp > 1), we find 
good agreement with the virial expansions |18| (see inset 
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FIG. 3. The contact of a nearly homogeneous sample (about 
30% of the trapped atoms probed), versus T/Tf at unitar- 
ity (black circles). The shaded area marks the superfluid 
phase transition, with some uncertainty in its exact position 
(T C /T F = 0.16 - 0.23) [25]. As a comparison, we plot the 
gaussian pair-fluctuation NSR model (GPF) [18], the self- 
consistent t-matrix model (GG) [25], the non-self-consistent 
t-matrix model (GoGo) [5], the 2nd and 3rd order virial ex- 
pansion [IS] , a quantum Monte-Carlo calculation (QMC) [TS] , 
and the contact extracted from a thermodynamic measure- 
ment done at ENS [26] . The error bars represent one standard 
deviation. The inset shows the high temperature behavior of 
the contact, where we find good agreement with the virial 
expansion. 



of figure[3]). For T/Tp < 0.4, our data do not agree fully 
with any of the many-body theories. It is worth noting, 
however, that two of the theories (GPF and GoGo) pre- 
dict a higher value for the contact above the superfluid 
phase transition than below, which may be consistent 
with observed sharp decrease near T/Tp = 0.16. We 
note that the predicted T c /Tp has some uncertainty, as 
indicated by the shaded region in figure [3] The non-self- 
consistent t-matrix model (GoGo) predicts an enhance- 
ment of about 50% in the value of the contact around T c 
[5], which the data do not show. We also do not observe 
an increasing trend in the contact for T > T c , in contrast 
to a recent QMC simulation |19j . 

As can be seen from Eq.Q, the contact is naturally 
normalized by Nkp, and the detuning by the Fermi en- 
ergy. However, a question which arises is how to de- 
fine Ep in our experiment. For a harmonically trapped 
gas, Ep is defined in terms of the trap parameters 
Entrap = MK^z) 1/3 (6A0 1/3 - On the other hand, the 
Fermi energy of a homogeneous gas is given in terms of 
its density (in one spin state), n: Ep^ om = |^-(67r 2 n) 2 / 3 . 
In our experiment, as we increase the power of the 
hollow light beams, we probe a smaller portion of the 
gas that is more homogeneous. The relevant Fermi 
energy, which we use in figures [2j and |3j is therefore 
the average of the local (homogeneous) Fermi energy: 
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FIG. 4. Contact versus the fraction of atoms probed for a gas 
with T/Tf = 0.46 at the center of the cloud. In the main 
plot, the measured contact (squares) is normalized in respect 
to the trap Icf, and is compared to the predictions of several 
theoretical models (lines) using the local density approxima- 
tion. The measured contact increases as we probe fewer atoms 
at the cloud center, where the local density is largest. The 
inset shows the contact normalized by the average Uf of the 
probed atoms (squares), compared to theoretical predictions 
of the homogeneous contact at the average T/Tf (lines). 



E 



_F,avg 



2mN, 



J P(r)n(r)[67r 2 n(r)] 2 / 3 d 3 r, where P(r) 
is the detection probability after optical pumping, and 
N p — J P(r)n(r)d 3 r is the number of atoms probed. 

The average local Fermi energy can be obtained from 
the density distribution of the atoms, n(r), and the de- 
tection probability, P(r). We measure n(r) by turning 
the trap off, without applying the hollow light beams, 
and imaging the cloud after 4 ms of expansion at the res- 
onance. To determine the density distribution in trap, we 
fit the distribution measured after expansion and rescale 
the dimensions back to t = 0, assuming hydrodynamic 
expansion [211 [27]. For the fit, we use the Thomas- Fermi 
distribution, which is known to fit the data well [2J. 

We obtain P(r) using a model of the optical pumping 
by the hollow light beams [22] . In the model, we assume 
that atoms that scatter a single photon are transferred to 
the dark state, and we account for the attenuation of the 
hollow light beams as they propagate through the cloud. 
For a given n(r), the propagation model gives us P(r) 
after the consecutive application of the two hollow light 
beams. We note that the results presented in figure [3] are 
not sensitive to the details of the model [24"] . 

In figure [4j we show the contact at T/Tp = 0.46 as a 
function of the fraction of atoms probed, which is var- 
ied by changing the intensity of the hollow light beams. 
The main part of figure [4] shows the contact per par- 
ticle in units of fc^trap m order to show the change 
in the measured signal. We find that the signal in- 
creases as we probe fewer atoms near the center of the 
trapped gas. We compare our results with several theo- 
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retical models, where the model lines are calculated by 
J P(r)n(r)C^[T/T F (v)}k F (v)dh, 
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(T/Tp) being the model prediction for a 
homogeneous contact (normalized to Nkp), Tp(r) — 
Ep(v)/kB is the local Fermi temperature, and fee is the 
Boltzmann constant. We find good agreement of the data 
with the models. 

In the inset of figure |4j we plot the contact di- 
vided by the average local kp, defined in the same 
way as in figure [3j For comparison, we also plot the- 
ory predictions for the homogeneous contact at the av- 
erage T/Tp, C^ cl {(T/T F )), where the notation () 
stands for density-weighted averaging. A reasonable 
criterion for homogeneity is when C™°^ cl ((T/Tp)) 
(C^ cl (T/7>)). When the fraction of the atoms probed 
is less than 30% we find that this approximation holds 
to better than 2% [2J]. When probing 30% of the atoms, 
we calculate that the rms spread in the local Tp has 
been reduced to about 20%. We find that the data for 
T/Tp = 0.46 and fractions lower than 30% agree with 
theory predictions for a homogeneous gas (see inset of 
figure [4]) . 

Lastly, we describe our determination of the tempera- 
ture of the gas at unitarity. Thermometry of a strongly 
interacting gas is not trivial, and different groups have 
used various techniques, including thermometry with a 
minority component 26 , measurement of the energy ver- 
sus entropy relation [28], and an empirical temperature 
extracted from fitting the cloud to a Thomas-Fermi dis- 
tribution |20j . We base our thermometry on a measure- 
ment of the release energy of the gas and the recently 
reported equation of state [29] . We determine the re- 
lease energy by taking an image of cloud after 4 ms of 
expansion at unitarity. Knowing our trapping potential, 
the equation of state, and the generalized virial theorem 
at unitarity [28] . we are left with only the temperature, 
T, as a free parameter in the calculation of the release 
energy. We find T by matching the calculated energy to 
the measured one [24]. We estimate that the one sigma 
uncertainty in the temperature is 5%. When reporting 
T/Tp in figure[3j we use Tp — Ep . avg /fcB- 

In summary, we have presented a measurement of the 
homogeneous contact of a unitary Fermi gas versus tem- 
perature. Our measurement is based on a novel technique 
that allows us to probe local properties of the cloud. Our 
data show good agreement with theory predictions for 
T/Tp > 0.4, but at lower temperatures no single predic- 
tion fully agrees with the data. Furthermore, the data do 
not show an enhanced narrow peak around T c , which was 
predicted to exist due to pair fluctuation in a pseudogap 
phase. To provide additional insight into the nature of 
the normal state of the unitary Fermi gas, it will be inter- 
esting to test directly the pseudogap pairing instability 
by combining our probing technique with momentum- 
resolved rf spectroscopy [2TJ [23] . 
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I. SETTING THE RF POWER FOR THE LINE SHAPE MEASUREMENTS 



The rf line shape is taken such that the number of atoms outcoupled by the rf pulse is 
small compared to total number of atoms. In order to determine the rf power that complies 
with this requirement, we have measured the number of atoms transferred by the rf pulse 
as a function of the rf power for different frequencies. The data presented here was taken 
without the optical pumping beams, for a gas at unitarity with ~ 80000 atoms per spin 
state and Ef/h = 7700 Hz, with h being the Planck constant. The temperature of the gas 
T/Tp = 0.13 was measured before ramping to the resonance field. 




(rf power) 1 ' 2 [arb] rf power [arb] 

FIG. SI. (a) Coherent oscillations with increasing rf power at a detuning of 5 MHz relative to the 
single-particle transition frequency. To show the coherent oscillations, the x-axis is the square root 
of the rf power (which is proportional to the Rabi frequency), (b) Saturation with increasing rf 
power at a detuning of 5 MHz (zoomed in on the initial rise in figure SI, plotted as function of the 
rf power). The solid (red) line is the fit to the exponential model introduced in the text. 

For small detunings, we observe oscillatory behavior reminiscent of coherent Rabi oscil- 
lations (figure Sla). For large detunings, the transfer is incoherent (see figure S2), and 
the number of outcoupled atoms saturates to a value that is between a quarter and a 
half of the total number of atoms (depending on the rf detuning). We fit the data to 
iV ou t = A(l — e~ p / p °) to find the saturation power Pq, and the results are given in Table I. 
For the lowest detuning, where coherent effects exist, we fit the signal up to the first peak 
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FIG. S2. Saturation with increasing rf power at a detuning of 45 MHz. The solid (red) line is the 
fit to the exponential model introduced in the text. 

(figure Sib). When measuring the rf line shape we keep the power less than Po/5, where the 
number of atoms transferred by the rf pulse is small compared to the total number of atoms. 
In this range, AUt ~ A^- to within 10%. We use a larger rf power at larger detunings and 
then linearly scale the number of atoms at each detuning to the power used for the central 
part of the line shape. 



detuning [kHz] 


P [arb] 


5 


0.33 


25 


2.5 


45 


8.4 


65 


12.7 


85 


17.1 



TABLE I. The measured saturation power, Pq, in arbitrary units, for different rf frequency detun- 
ings. 
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II. OBTAINING THE IN-SITU DENSITY DISTRIBUTION 



We use the in-situ density distribution, n(r), in order to calculate the average kp and 
E F of the probed atoms. To get n(r), we use the fact that at unitarity the cloud expands 
hydrodynamically, and the dynamics are governed by the continuity equation. The solution 
for the continuity equation with harmonic confinement with a time-dependent trapping 
frequency ou(t) is self-similar with the following scaling transformation: rj(t) = 6j(t)rj(0), 
where r\ is the spatial coordinate [i = x,y,z), and bi(t) obeys the equation [1, 2]: 

S.W = -Mt> + m]b ^ M • ID 

with the initial conditions 6^(0) = 1 and 6«(0) = 0. The constant 7 is the characteristic 
exponent in the equation of state /x(n) oc n 7 , where /z is the chemical potential and 7 = 2/3 
at unitarity. For a sudden turn off of the trap, Wj(0) is the trapping frequency along the % 
axis, and uii{t > 0) = 0. 

In figure S3, we plot the measured width of the cloud in the axial and radial directions 
as a function of the expansion time. The data show a rapid increase in the size in the radial 
(tight) direction of the cloud and almost no increase in the axial direction-a characteristic 
of hydrodynamic expansion. The solid lines show the numerical solution of Eq.(l) with 
oj r = 2n x 226 Hz and oj z = 2ir x 19 Hz (which were measured independently), which 
agrees very well with the data. We find that after 4 ms of expansion the finite resolution 
of the optical system does not affect the extracted parameters, and therefore we choose this 
expansion time for the density measurements. We fit the measured density profiles at 4 ms 
with a Thomas-Fermi distribution, which we find to be general enough for this purpose. 

We have tested our density determination method by looking at the ratio of the peak 
density at unitarity to the peak density of a weakly interacting gas, at low temperatures. 
The density distribution of the weakly interacting gas is measured using a fit to a Thomas- 
Fermi distribution after ballistic expansion. For T = 0, this ratio is ny /no = £ -3 / 4 , where £ 
at unitarity is a universal constant that relates the chemical potential to the Fermi energy: 
H = ^ep. From the measured density ratio at a temperature of T/Tp trap = 0.15, we extract 
a value of £ = 0.40 ± 0.05, which is consistent with recent determinations of this universal 
constant by other groups [3, 4]. 
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FIG. S3. Hydrodynamic expansion at unitarity. We start with a weakly interacting gas with 
~ 90000 atoms per spin state at T/Tp = 0.12 and ramp adiabatically to the Feshbach resonance 
field. We fit the cloud with a Thomas-Fermi distribution after a variable expansion time and 
extract the rms widths, cr(t), in the radial and axial directions. For the data, we de-convolve the 
measured width with a gaussian point spread function with an rms width of 2.9 fim, to account 
for the finite resolution of the optical system. The data is normalized by the initial cloud size, 
which is 33.4 /im and 2.8 /im in the axial and radial directions, respectively. The solid lines are 
the numerical solution of the hydrodynamic equation. 

III. THERMOMETRY 

Our thermometry assumes a knowledge of the trapping potential V(r) and the equation 
of state n(/i, T), where T is the temperature and \x is the chemical potential. For a non- 
interacting gas the equation of state is known, and for the unitary gas, we use the equation 
of state recently measured at MIT [3]. The trapping potential is calibrated from the known 
optical trap beam waists and the measured oscillation frequencies in all three directions. 
We adopt a local density approximation approach; the local chemical potential is given by 
/i(r) = fiQ — V(r). For a given T and number of atoms, N, /io is set by the normalization 
requirement N = J n\p(r), T] d 3 r. The equation of state then determines the complete 
density profile n(r), from which we can calculate other quantities such as the entropy, total 
energy, release energy, and shape of the cloud. Since we measure N, we get a one-to- 
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one correspondence between T and these quantities, and therefore any of them can serve 
as a thermometer. With the unitary gas, we have chosen to use the release energy as a 
thermometer. 

We determine the release energy by measuring the cloud after it expands at unitary for 
4 ms. The release energy per particle is calculated from the measured density profile of the 
expanded gas using 

E rel = E ^ = jf E /f (y) 2 K«-n (r)]d 3 r , (2) 

i=x,y,z i=x,y,z 

where is the corresponding spatial coordinate (i = x,y,z), t is the expansion time, and 
n t (r) is the density distribution at time t. In the experiment, we use t — 4 ms. We have 
verified that the release energy measured at t = 4 ms is the same as that measured after 12 
ms of expansion. For a given potential V(r), the release energy is given by [5]: 

£ rc i = ^(r-W(r)> , (3) 

where the symbol (} stands for the density-weighted average: (g(r)) = jf J g(r)n(r)d 3 r. By 
equating the calculated E re i(T) to the measured E m \, we determine T. 

As a comparison, we have used two other techniques to extract the temperature of the 
unitary gas. The first technique we compare to is based on the widely used practice of 
fitting the strongly interacting gas to a Thomas-Fermi distribution and extracting an em- 
pirical temperature, T, from the fitted fugacity [6]. At T = 0, the empirical temperature is 
connected to the real temperature by T = Ty^, where £ is the universal constant defined 
above [2]. Albeit without a complete theoretical justification, one can then extend this to 
finite temperatures and extract T [6]. In the following analysis, shown as the blue triangles 
in figure S4, we used £ = 0.376 [3]. 

The second thermometry method we compare to is based on the entropy of the weakly 
interacting gas before the ramp to unitarity. We calculate the entropy of the weakly interact- 
ing gas from the measured temperature and the trapping potential. In the experiment, we 
start from the weakly interacting gas and slowly ramp to the Feshbach resonance field. By 
performing this ramp there and back and comparing the entropy before and after the ramp 
for a gas initially at T/T F = 0.12 and T/T F = 0.22, we have determined that the entropy 
increases by about 6% when going to the Feshbach resonance field. Assuming this increase, 
we use the entropy of the unitary gas together with the equation of state as the thermometer. 
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In figure S4, the T/T F we obtain from these two additional techniques are plotted against 
the release energy thermometry. We find a good agreement between all the three techniques 
up to T/Tp = 0.4. Above that temperature, the empirical temperature technique becomes 
unreliable since the effect of quantum degeneracy of the shape of the cloud diminishes. The 
entropy technique starts to show a small systematic deviation upwards above T/Tp = 0.4. 
The close agreement of the three techniques, which are based on independent observables, 
up to T/Tp = 0.4 gives us confidence in our thermometry. To estimate the errors in T we 
look at the difference between the entropy and release energy thermometry techniques. 
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FIG. S4. Comparison of different thermometry methods. The x-axis is the temperature, T/Tp, 
where Tp is the trap Fermi temperature, as extracted from the release energy. The y-axis is the 
temperature we get from the two other thermometry methods (see text for more details). The 
dashed line is y = x. 



IV. THE EFFECT OF THE REMAINING DENSITY INHOMOGENEITY 

In figure 3 in the paper, we compare our data to the predictions of several theoretical 
models. Here we show that the effect of the remaining density inhomogeneity of the probed 
sample on the theory predictions is negligible. We use the detection probability, P(r), and 
the density distribution, n(r), for each of the data points, to calculate the average contact 
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predicted by each theoretical model according to: 

(C) = J P(r)n(r)C^ cl [T/T F (r)]k F (r)dh , (4) 

where C™^ el is the prediction for the contact of a homogeneous gas theory (normalized to 
Nkf), N p is the number of probed atoms, and (kp) is the average k F . The comparison of 
the average contact, (C), and the homogeneous contact for three different models is shown 
in figure S5. The graph clearly demonstrates that (with 30% of the atoms probed) the effect 
of the remaining density inhomogeneity on the contact is negligible, and hence theories for 
the homogeneous contact can be compared directly to the data. 
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FIG. S5. Comparison of the homogeneous contact calculated by several theoretical models (solid 
lines) and the contact averaged over the remaining density inhomogeneity when probing the central 
30% of the cloud for the same models (open symbols). The excellent agreement of the points and 
the lines shows that the effect of the remaining density inhomogeneity on the contact data can be 
neglected. 

V. THE OPTICAL PUMPING MODEL 

To calculate the spatially dependent probability that atoms are probed, -P(r), we use a 
model for the optical pumping by the hollow light beams. In Ref. [7] we have introduced 
such a model, which was also tested experimentally The calculated -P(r) along with n(r) is 
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used in our determination of the average quantities (density, k F , T F , E F and the contact) of 
the probed atoms. We show here that the results presented in this paper are not sensitive 
to the details of the optical pumping model. 



4.0 



3.5 - 




■ full removal model 
* sharp cut model 



3.0 - 



* 2.5 - 
O 2.0- 



1.5 - 



1.0 



0.0 



0.5 



1.0 




FIG. S6. The homogeneous contact at unitarity as function of the temperature, using the full 
removal model (squares) and the sharp cut model (triangles). The good agreement of the two 
models shows that the exact shape of the detection probability, -P(r), does not affect the contact 
data. This is because the central part of the trapped gas is nearly homogeneous, and the data was 
taken probing the central 30% of the atoms. 

To show this, we introduce a second, much more simplified model, which we refer to as 
the "sharp cut" model. In the sharp cut model, we assume that atoms are left for probing 
only if their position (x, y, z) satisfy x 2 + y 2 < R\ and z 2 + y 2 < R 2 , where R\ and R 2 are set 
to reproduce the fraction of atoms probed after the application of one or both of the hollow 
beams. In figure S6, we present the homogeneous contact, similar to figure 3 in the paper, 
but using the two different models. As can be clearly seen in the figure, the results of the two 
models are essentially the same, which demonstrates that the details of the calculation of 
P(r) are not important when the fraction of atoms probed is as small as 30%. In figure S7, 
we show the average Fermi energy for the data, where again the two models give consistent 
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FIG. S7. The average Fermi energy with ~ 30% of the atoms probed, using the full removal model 
(squares) and the sharp cut model (triangles). 
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